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ABSTRACT 

Sixteen pulsars have been discovered so far in blind searches of photons collected with the Large Area 
Telescope on the Fermi Gamma-ray Space Telescope. We here report the discovery of radio pulsations from 
two of them. PSR J1741-2054, with period P = 413 ms, was detected in archival Parkes telescope data and 
subsequently has been detected at the Green Bank Telescope (GBT). Its received flux varies greatly due to 
interstellar scintillation and it has a very small dispersion measure of DM = 4.7 pc cm" 3 , implying a distance 
of k 0.4 kpc and possibly the smallest luminosity of any known radio pulsar. At this distance, for isotropic 
emission, its gamma-ray luminosity above 0.1 GeV corresponds to 25% of the spin-down luminosity of E = 
9.4 x 10 33 ergs _1 . The gamma-ray profile occupies 1/3 of pulse phase and has three closely-spaced peaks 
with the first peak lagging the radio pulse by <5 = 0.29 P. We have also identified a soft Swift source that is 
the likely X-ray counterpart. In many respects PSR J1741-2054 resembles the Geminga pulsar. The second 
source, PSR J2032+4127, was detected at the GBT. It has P = 143 ms, and its DM = 115pccm" 3 suggests 
a distance of « 3.6 kpc, but we consider it likely that it is located within the Cyg OB2 stellar association at 
half that distance. The radio emission is nearly 100% linearly polarized, and the main radio peak precedes by 
S = 0.15P the first of two narrow gamma-ray peaks that are separated by A = 0.50P. The second peak has 
a harder spectrum than the first one, following a trend observed in young gamma-ray pulsars. Faint, diffuse 
X-ray emission in a Chandra image is possibly its pulsar wind nebula. PSR J2032+4127 likely accounts for the 
EGRET source 3EG J2033+4118, while its pulsar wind is responsible for the formerly unidentified HEGRA 
source TeV J2032+4130. PSR J2032+4127 is coincident in projection with MT91 213, a Be star in Cyg OB2, 
although apparently not a binary companion of it. 

Subject headings: gamma rays: observations — ISM: individual (TeV J2032+4130) — open clusters and asso- 
ciations: individual (Cyg OB2) — pulsars: individual (PSR J1741-2054, PSR J2032+4127) 
— X-rays: individual (Swift J174157.6-205411) 
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1. INTRODUCTION 

Neutron stars hold vast stores of rotational energy. Mag- 
netic braking in such stars generates spin-down luminosities 
of E « 10 28 " 38 ergs" 1 . A minute portion of this emerges in 
collimated coherent radio beams through which most known 
rotation-powered pulsars are detected. Much in pulsar elec- 
trodynamics remains obscure, including processes and loca- 
tions of particle acceleration in the magnetosphere. In con- 
trast to radio emission, pulsed gamma-ray luminosity can be a 
substantial fraction of E, and its study therefore holds promise 
for significant advance in understanding magnetized rotating 
neutron stars. 

Until recently, only six p ulsars were conf irmed gamma- 
ray emitters above 0.1 GeV dThom pson 2004). Their detec- 
tion in the 1990s with the EGRET experiment on the Comp- 
ton Gamma-Ray Observatory used rotational ephemerides to 
fold the sparse gamma-ray photons at a known period. The 
ephemerides were obtained from radio observations, with 
the exception of Geminga, which is not a detected radio 
source. Given these limitations, it was not possible, for in- 
stance, to establish the fraction of gamma-ray pulsars that do 
not emit radio beams detectable at the Earth, i.e., that are 
"radio quiet". Acc ording to predictions of oute r magneto- 
sphere models (e.g., iRomani & Yad igaroglu 19951) . gamma- 
ray beams are broader than radio beams and thus are poten- 
tially detectable from a larger solid angle. Measuring the frac- 
tion of radio-quiet gamma-ray pulsars therefore can provide 
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FIG. 1.— Radio detections of PSR J1741-2054: Parkes at 1.4 GHz (left); 
GBT at 2 GHz (center), and at 350 MHz (right). In each panel the pulse 
profile is shown repeated in phase, as a function of time in the grey scale, and 
summed at the top. The relative pulse phases between panels are arbitrary. 

information on beam shapes and emission regions, an d con- 
strain emission mechanisms (e.g jHarding et all2007l) . 

The Fermi Gamma-ray Space Telescope, operational since 
mid-2008, is being used to make a torrent of pulsar discover- 
ies with its Large Area Telescope (LAT), which has substan- 
tially greater effective area, field of view, and angular resolu- 
tion than EGRET. Using radio-derived ephemerides, at least 
14 pulsars have been detected for th e first time by LAT , of 
which eight are millisecond pulsars dAbdo et alj|2009b1lcllg l). 
In addition, at least 16 pulsars have been discovered in peri - 
odicity searches of gamma-ray photons dAbdo et al.l l2009f ). 
While a few of these had been searched previously in very 
sensitive rad io observations, and therefore may be considere d 
radio quiet (iBecker et all 12004 lHalpern et all l2004 l2007t) . 
most new LAT pulsars should be searched for radio counter- 
parts before their true character can be determined. A radio 
detection yields an estimate of distance and therefore of lu- 
minosity, allowing for a more comprehensive study of these 
energetic neutron stars. Here we report the detection of radio 
pulsations, as well as additional gamma-ray and X-ray analy- 
sis, from two Fermi LAT pulsars. 

2. OBSERVATIONS AND RESULTS 

Nine of the new LAT pulsars are located in areas pre- 
viously surveyed with the ATNF Parkes radio telescope. 
Three had been targe ted in sensitive but unsuccessful searches 
dRoberts et alJ |2002|) . We searched anew these data, as 
well as archival 1.4 GHz data for the other six LAT pul- 
sars (J0633+0632, J 1 459-60, J 1 732-3 1 J 1 74 1-2054, J 1 8 1 3- 
1246, and J1907+06: lAbdo et alJl2009l) . The 95% C.L. er- 
ror radius for these LAT sources spanned 4'6— 8 '4. With a 
beam radius of 7', there were in general multiple relevant 
Parkes pointings to inspect. We searched eac h such data set 
in dispersion measu re (DM) with PRESTO ( Ransom] EMU 
iRansom et al.l20 02). folding time samples from the 35 minute 
integrations modulo the pulsar period predicted from each 
LAT ephemeris. In the case of PSR J1741-2054, we detected 
pulsations. 

2.1. PSR J 1741-2054 
2.1.1. Radio measurements 



The 1.4GHz detection of the P = 413ms PSR J1741- 
2054 (left panel of Figure Q} is from data collected across 
a bandwidth of 288 MHz in the Parkes multibea m survey 
of the Galactic plane (e.g., iManchester et all 1200 lb on 2000 
November 24. This location (RA. = 17 h 41 m 5L s 3, deck = 
-21°01'10"; all positions here refer to the J2000.0 equinox) 
was 1'.6 away from the nominal LA T position, which had an 
uncertainty of 84 (lAbdo et alj|2009d) . 

In order to further constrain its position, we attempted sev- 
eral observations of the pulsar with the NRAO Green Bank 
Telescope (GBT) at 2 GHz, for which the beam radius is 3', 
using GUPPI 21 with approximately 700 MHz of usable band- 
width. We did this at six epochs, each time on a grid of point- 
ings surrounding the nominal Parkes position. We only de- 
tected the pulsar on 2009 March 7, at RA. = 17 h 41 m 50?l, 
decl. = -20°55'02", 4.'6 from the LAT position. The received 
flux varied greatly during the 5 minute observation (middle 
panel of Figure [TJ, possibly due to interstellar scintillation. 

We obtained from the first two detections only an up- 
per limit on DM. We finally measured the dispersion of 
PSR J1741-2054 with a 350 MHz detection on March 12 
(right panel of Figure [B using the GBT BCPM dBacker et alJ 
119971) with a bandwidth of 48 MHz: DM ^5pccm" 3 . Such 
a low DM points to scintillation (see, e.g.. lLorimer & Kramer! 
l2005h as an explanation for the several non-detections (which 
further include at Parkes two observations at 1 .4 GHz and one 
at 3 GHz, and at GBT one 0.8 GHz observ ation). For this 
DM an d location at (l,b) = (6.°42,+4.°90), the lCordes & Laziol 
(2002) electron density model yields a pulsar distance of 
d = 0.4 kpc, with significant uncertainty. Hereafter, we pa- 
rameterize this distance by c/o.4 = c//(0.4kpc). 

We have estimated the observed flux densities S„ for the 
three radio detections: So.35 ~ 1.33 mJy, Si a « 0.16mJy, and 
S2 ~ 0.09 mJy. However, because the detection positions at 
the two highest frequencies are slightly offset from the true 
position (e.g., Figure |2), these values somewhat underesti- 
mate the flux on those days. Conversely, because of scin- 
tillation, the average flux is likely substantially lower than the 
position-corrected values on those days. Thus we cannot infer 
a spectral index. 

A comparison of the barycentric pulsar periods determined 
from radio observations taken more than 8 years apart implies 
a period derivative of P = (2.7 ± 1.1) X 10~ 14 (unless otherwise 
noted, all uncertainties here are given at the 1 a C.L.). 

2.1.2. Gamma-ray and radio timing and position 

Using LAT data spanning one year, we have extracted 17 
gamma-ray times-of-arrival (TOAs) for PSR J1741-2054 that 
we use to obtain a timing solution with TEMPO 22 . While the 
fitted position in R.A. has an uncertainty of only 1", that in 
decl. is uncertain by 0.'6 (see Figure 0, owing to the low 
ecliptic latitude of the pulsar. In addition, a Swift observation 
reveals a probable X-ray counterpart at R.A. = 17 h 41 m 57. s 6, 
decl. = -20°54'11", with uncertainty of about 4". The Swift 
source position is 5" in R.A. and 58" in decl. from the tim- 
ing position, consistent within the combined 2 a uncertainties. 
Because of the still significant positional uncertainty (which 
should be essentially eliminated with a planned Chandra ob- 
servation), we cannot rule out the presence of a small amount 
of timing noise. In any case, the timing solution already yields 
an accurate P = 1 .69 x 10" 14 . In turn, t c = P/(2P) = 0.39 Myr 

21 https://wikio.nrao.edu/bin/view/CICADA/GUPPiUsersGuide 

22 http://www.atnf.csiro.au/research/pulsar/tempo 
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FIG. 2.— X-ray field of PSR J1741-2054: portion of Swift 0.2-10keV 
image smoothed with a a = 7" Gaussian. The soft source Swift J174157.6- 
20541 1, with 69 photons within an extraction radius of Of 5, is the only plau- 
sible counterpart in this 4.4ks observation for PSR J1741-2054, whose 3<r 
timing position (ellipse), 95% C.L. gamma-ray point source (solid circle), 
an d FWH M radio (dashed circle) positional uncertainties are indicated (see 
jjj) |2T2l and |2T4l . 



and£ = 9.4x 10 33 ergs _1 . 

We determined the position for the PSR J1741-2054 LAT 
point source with the binned likelihood application ptlike, 
by selecting photons above 0.125GeV from the light curve 
peak, with phases 0.24 < </> < 0.55 (Figure [3j this phase 
range subtends roughly the half-maximum level of pulsed 
emission, offering the best signal-to-noise ratio for localiza- 
tion), and varying the position to maximize the point source 
significance. The resulting position is R.A. = 17 h 41 m 53 ? 1 , 
decl. =-20°54'51", with a 95% C.L. error radius of 2(3. This 
is only l'.2 away from the Swift source. The best radio posi- 
tion is also consistent with the LAT timing and Swift source 
positions (see Figure [2] for an illustration of the various posi- 
tional determinations). 

Using the 350 MHz and 2 GHz TOAs in a separate fit, 
we also obtain a higher-precision measurement of disper- 
sion, DM = (4.7 ±0.1)pccm" 3 . This allows the determina- 
tion of the phase offset between the radio and gamma-ray 
pulses separately for both recent radio detections. The re- 
sults are identical, and we show in Figure |3]the phase relation 
for the 350 MHz and gamma-ray profiles. Th e LAT profiles 
are co nstructed from "diffuse class" photons (lAtwood et al.l 
2009) selected from a region with an energy-dependent ra- 
dius OPSOE/lGeV) -0 ' 75 (with maximum and minimum radii 
of 0.°8 and 0.°35, respectively) around the best position over 
the time range 2008 August 4 to 2009 July 4. We used ener- 
gies from 0.2-10 GeV. The phase alignment was done using 
the gtpphase 23 tool to assign phases to each LAT photon ac- 
cording to the best timing model. We found that the pulsed 
emission was well-described by a three-Gaussian model and 
extracted the parameters with a maximum likelihood fit to 
the unbinned profile photon phases. We applied the model 
to the energy bands shown in Figure [3] (0.2-1.0 GeV and 1- 
lOGeV) as well as to the full (0.2-10GeV) energy range. We 

23 http://fermi.gsfc.nasa.gov/ssc/data/analysis 
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Pulse Phase 

FIG. 3. — Phase-aligned Fermi and GBT pulse profiles of PSR J1741— 
2054. In the upper panel, the radio profile is displayed with an arbitrary 
intensity scale along with LAT counts in the 0.2-1.0 GeV band. The bottom 
panel shows the higher-energy LAT counts and a comparison between the 
two panels shows clear evolution of the peak structure (PI weakening, P3 
strengthening) with energy. At the highest energies (> 3 GeV), P3 dominates. 
The displayed gamma-ray and radio profiles have, respectively, 32 and 128 
bins per period. Two full rotations are shown. 



found no significant difference in the widths and positions 
of the Gaussians, so we report here the results from a fit to 
the full energy range. The first gamma-ray peak ("PI"), with 
phase 0.292 ± 0.005 and FWHM = 0.06 ± 0.01, determines 
the phase offset with the radio peak, S = 0.29 ± 0.02. The sec- 
ond gamma-ray peak ("P2") occurs at <f> = 0.384 ± 0.015 with 
FWHM = 0. 1 3 ± 0.03, while P3 is located at = 0.5 1 8 ± 0.015 
with FWHM = 0.11 ±0.02. Interpreting the P3-P1 separa- 
tion as the conventional separation in a two-peak profile yields 
A = 0.226 ±0.016, although it is unclear whether this is mean- 
ingful for such a complex profile. 

2.1.3. Gamma-ray spectrum 

To fit the gamma-ray spectrum of PSR J1741-2054, we 
have performed an unbinned likelihood analysis with the ap- 
plication gtlike 23 , and an additional check usin g the inde- 
pendent ptlike. Both approaches are outlined in lAbdo et alj 
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(2009d), and here we use an updated instrument response 
function, P6_v3_dif f 23 , that corrects a pileup effect iden- 
tified in orbit. The on-axis energy resolution of LAT is 
(equivalent Gaussian la) 15%— 9% in the 0.1-1 GeV range, 
8% over 1-lOGeV, and 8%-18% in the 10-300GeV range 
(lAtwood et al.l 12009). We used gll_iem_v0 2 23 as the 
model for the Galactic diffuse background. We selected 
gamma-rays with energy > 0.1 GeV and with zenith an- 
gle < 105°, extracted from within a radius of 15° around 
the pulsar over the time range 2008 August 4 to 2009 July 
4. We modeled all significant point sources in the aper- 
ture with a power-law spectrum. We found no evidence of 
point source emission off pulse (0.68 < <f> < 1.18), so to 
increase the signal-to-noise ratio, we considered only pho- 
tons with phase 0.18 < <f> < 0.68. The likelihood ratio 
test indicates that an exponentially cut-off power law given 
by dN/dE = N (E/l GeV) _r exp(-E/E c ) is preferred over a 
simple power-law model with 12(7 significance. We find a 
photon index V = 1.4 ±0.1 ±0.1 and a cut-off energy E c = 
(1.1 ±0.2 ±0.2) GeV, giving a photon flux above 0.1 GeV of 
(20 ± 1 ± 3) x 10 _8 cm- 2 s" 1 and an energy flux of F 7 = (12 ± 
1 ± 2) x 1 _1 1 erg cm" 2 s -1 , where the first (second) errors in- 
dicated are statistical (systematic). The resulting gamma-ray 
luminosity, assuming effectively isotropic emission (/h =1), 
is L 1 = AnfaF^d 2 = 2.4 x lO^c/^ergs -1 = 0.25 dl A E. 

The third peak has a significantly harder spectrum than 
PI: it is not detectable at 0.1-0.3 GeV, then becomes grad- 
ually more significant with increasing energy, such that above 
1 GeV its peak count rate exceeds that of PI. The spectrum 
of P2 appears to be of intermediate hardness between that of 
P3 and PI (see Figure |3}. The highest energy photon likely 
detected so far from this pulsar arrived at phase 0.491 (in P3) 
with energy 5.8 GeV. 

2. 1 .4. Spectrum of likely X-ray counterpart 

The error circle of PSR J1741-2054 was observed with 
the Swift X-ray telescope on 2008 October 16, obtaining 
4.4 ks of exposure with 2.5 s time sampling. The soft source 
Swift J 174 157. 6-2054 1 1 near the field center is unresolved 
and provides 69 background-subtracted counts in the 0.3- 
5keV range. There is no other detected X-ray source that 
is compatible with the position of the pulsar (Figure |2), and 
the spectral character of this source (see below) is compatible 
with our expectations for PSR J1741-2054. We therefore be- 
lieve that Swift J174157.6-205411 is the likely pulsar coun- 
terpart. 

We see no evidence in Swift J174157. 6-205411 for inter- 
stellar absorption. If we fit the source spectrum with a fixed 
Na = 1.5 X 10 20 cm" 2 (10 times the DM-determined free elec- 
tron column), we obtain a power-law index T = 2.5 ± 0.4 
and F x (0.5- lOkeV) = 6.4 x lO^ergcrrrV 1 . Allowing N H 
to vary yields Nh < 3 x 10 21 cm" 2 and a poorly determined 
T = 2.3 ±1.7. Thermal models give kT = (0.2 ±0.1)keV but 
are not well constrained because of the unknown absorption. 

The observed X-ray flux likely represents a composite spec- 
trum similar to other relatively old and nearby gamma-ray 
pulsars such as Geminga and PSR B 1055-52. At an age com- 
parable to the characteristic age of PSR J1741-2054, 0.4 Myr, 
thermal emission from the full surface is expected to have a 
relatively low kT « 50-70eV. Such low temperatures alone 
do not reproduce the Swift spectrum which has significant 
counts above 1 keV. Thus the observed X-rays likely have a 
significant synchrotron component coming from secondary 
e pairs that are gamma-ray annihilation products. 
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FIG. 4.— Polarimetric profiles of PSR J2032+4127 based on GBT observa- 
tions with GUPPI: 0.5 hr at 2 GHz (top) and 2.0 hr at 0.8 GHz (bottom). The 
black traces corresponds to total intensity, while the red and blue lines corre- 
spond, respectively, to linear and circular polarization. The position angles in 
each upper panel have been converted to the pulsar frame, using the measured 
RM = (+215± l)radnT 2 . 

2.2. PSRJ2032+4127 
2.2.1. Radio measurements 

On 2009 January 5 we used the GBT to search for radio 
pulsations from LAT PSR J2032+4127. We recorded data for 
1 hr from an 800 MHz band centered on 2 GHz, using GUPPI 
to sample each of 1024 frequency channels every 0.16 ms. 
The telescope beam, with radius 3', was centered on the LAT 
source at R.A = 20 h 32 n T3. s 9, d eck = +41°22'34" with error 
radius of 5.'1 (lAbdo et alJ]2009d) . 

We did a search in DM while folding the data modulo the 
71.6 ms period of the original LAT ephemeris, and detected 
the pulsar. Subsequently we did a blind search of the radio 
data and established that its true period is P = 143 ms. Fol- 
lowing some gridding observations, we obtained an improved 
position: R.A. = 20 h 32 m 14 s , decl. = +41°27'00", with uncer- 
tainty radius of 0.'7. At this accurate position, the pulsar is 
clearly detectable at the GBT in only 1 minute. 

We have also obtained two calibr ated polarimetric o bser- 
vations (analyzed with PSRCHIVE: iHotan et al.ll2004l) . one 
each at 2 GHz and 0.8 GHz, which show that the profile is 
nearly 100% linearly polarized, with a small interpulse de- 
tected at the lower frequency (Figure|4]i. The period-averaged 
flux density at 2 GHz is S2 = 0.12mJy, and Sq.% = 0.65 mJy. 
Each of these measurements has fractional uncertainty of 
about 20%, including an allowance for the small degree of 
interstellar scintillation inferred from many 2 GHz observa- 
tions. The resulting spectral index is a = -1.9 ±0.4, where 
S v cx v a . 
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FIG. 5.— Phase-aligned GBT and Fermi pulse profiles of PSR J2032+4127. 
The gamma-ray peaks are modeled as Gaussians of, respectively, 
FWHM/P = 0.026±0.003 and 0.051 ±0.005. The radio profile is displayed 
with an arbitrary intensity scale. The radio and gamma-ray profiles are dis- 
played with, respectively, 256 and 32 bins per period. Two full rotations are 
shown. 



coincides with the radio source and the optical star to within 
0"6, which is comparable to their combined statistical and 
systematic uncertainties. 

The timing fit also yields the phase offset between the ra- 
dio and gamma-ray profiles, which we show in Figure[5] The 
LAT profiles are constructed from diffuse class photons se- 
lected from a region of radius 0.°8 around the best position 
over the time range 2008 August 4 to 2009 July 4. We used 
only energies above 0.3 GeV to help reduce background con- 
tribution from nearby sources in this crowded field. The phase 
alignment was done using gtpphase to assign phases to each 
LAT photon according to the best timing model. The phase 
spacing between the two gamma-ray peaks (A = 0.50 ± 0.01) 
and the offset between the radio and first (PI) LAT peak 
(8 = 0. 15 ±0.01) were determined by fitting the unbinned pho- 
ton phases to a two-Gaussian model. 

We obtained the LAT position for PSR J2032+4127 with pt- 
like by selecting photons above 0.125 GeV from 0.12 < <f> < 
0.20 and 0.60 < <fi < 0.72 (Figure 01, and varying the posi- 
tion to maximize the point source significance. The resulting 
R.A. = 20 h 32 m 15?8, decl. = +41°26'17", with a 95% C.L. er- 
ror radius of is only 1.'2 away from the pulsar timing 
position. 



2.2.2. Gamma-ray and radio timing and position 

We are timing PSR J2032+4127 at the GBT. The mea- 
sured DM= (114.8 ±0.1) pc cm" 3 , at (l,b) = (80.°22,+l.°03) 
corres ponds to d = 3.6 kpc, according to the ICordes & Laziol 
(2002) model. Hereafter, we parameterize this distance by 
c/3.6 = af/(3.6kpc). The average radio TOA uncertainty is 3 
times smaller than the corresponding gamma-ray value, but 
we have a radio data span of only 6 months compared to 1 
year for Fermi. As a consequence, the best overall timing 
solution is obtained from a joint fit. Some rotational insta- 
bility is detectable in PSR J2032+4127 as timing noise. This 
is parameterized in TEMPO by the (non-stationary) second 
derivative of the rotation frequency, (—1.7 ±0.3) x 10~ 21 s" 3 . 
The measured P = 2.00 x 10" 14 implies r c = 0. 1 1 Myr and E = 
2.7 x 10 35 ergs _1 . The timing fit gives positional uncertain- 
ties in R.A. and decl. of 0."4 and 0"2, respectively, although 
timing noise may contribute a systematic error amounting to 
-1". 

Radio imaging of th is region was c onducted by 
iParedes et all (2007) and iMarti et alJ (12007b . including a 
6 10 MHz survey with the Giant Metrewave Radio Tele- 
scope (GMRT), to identify potential counterparts for the 
yery-high-energy gamma-ray source TeV J2032 +4130 
dAharonian et alJ l2002t 120051) . IParedes et alJ d2007l) noted 
that their GMRT sou r ce 5 c oincides with star 213 in the 
Mass ey & Thompson! (1199 ll hereafter MT91) survey of 
massive stars in Cyg OB2, as well as with a Chandra X-ray 
source. The latter associations are explored further in § 12.2.41 
Here we conclude that PSR J2032+4127 and GMRT source 
5 are one and the same based on positional coincidence (see 
Table [B, and because the (0.5±0.1)mJy flux of GMRT 
source 5 at 6 10 MHz is compatible within the uncertainties 
with the radio spectrum of PSR J2032+4127 (§ |2~2~TI Thus, 
the position of PSR J2032+4127 is known to 0."5 in each 
coordinate irrespective of timing noise. We are able to limit 
possible systematic uncertainties in the astrometry of the 
optical reference frame relative to the radio to w 0"1 using 
the Tycho position of the star Cyg OB2 4, which agrees with 
our optical astrometry to this level. Thus, the X-ray source 



2.2.3. Gamma-ray spectrum 

The Cygnus region contains several bright gamma-ray point 
sources and strong diffuse emission that remains difficult to 
model (see Figure [6}. However, as with PSR J1741-2054 
(§ 12.1.31 . we found no evidence for point source emission off 
pulse, and we found that we could measure the spectrum of 
PSR J2032+4127 robustly by selecting only photons from on- 
pulse phases (defined conservatively as being 0.08 <4>< 0.20 
and 0.56 < (j> < 0.80). We used the same time, energy, 
and zenith angle cuts as for PSR J1741-2054 (§ |2~T3T |. but 
with a 20° extraction radius, and modeled the bright nearby 
gamma-ray pulsars J202 1+3651 and J202 1+4026 with an ex- 
ponentially cut-off power-law spectrum. We find that an 
exponentially cut-off power law is preferred over a simple 
power law with 9 a significance. The best fit parameters, 
T = 1.1 ±0.2 ±0.2 and E c = (3.0 ± 0.6 ± 0.7) GeV, give a 
photon flux above 0. 1 GeV of (7 ± 1 ± 2) x 10" 8 cm -2 s" 1 and 
F 1 = (9 ± 1 ± 2) x 10 -11 erg cm" 2 s . For assumed isotropic 
emission, L 7 ~ 1.4 x 10 35 J 2 6 erg s" 1 = Q.5d\ b E. 

As with PSR J1741-2054 (§ I2J31 the trailing peak has 
a harder spectrum than PI: the P2/P1 ratio of baseline- 
subtracted maximum count rate is 1.0 at 0.3-1 GeV but rises 
to 1 .5 above 1 GeV, while the ratio of total counts in the peaks 
changes from 2.3 to 2.8 in the same energy ranges. The high- 
est energy photon likely detected so far from this pulsar ar- 
rived at phase 0.642 (in P2) with energy 9.8 GeV. 

2.2.4. X-ray and optical observations 

PSR J2032+4127 is located in the direction of the massive 
Cyg OB2 stellar association. Several X-ray studies have been 
made of the region containing the pulsar dButt et alJ 120061 : 
iHorns et all 120071: iMukherjee etail I2003L |2007j), to identify 
potential counterparts for TeV J2032+4130. Here we have 
reanalyzed a 50 ks Chandra X-ray Observatory archival ob- 
servation of the field containing PSR J2032+4127, obtained 
on 2004 July 12 with a ti me resolution of 3 . 2 s (observational 
detail s can be found in iButt et all l2006t [Mukherie e et alJ 
120071) . The pointing of this ACIS-I CCD observation placed 
PSR J2032+4127 3' off-axis, where the angular resolution is 
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FIG. 6. — The observed Fermi-LXI weighted counts map of the Cygnus region (20° X 20° centered on PSR J2032+4127). The photons, selected with the 
same time and zenith angle cuts used in the spectral analysis, and with energies between 0.2 and 10 GeV, were binned in energy and in position with a data pixel 
size appropriate to the point spread function (PSF) at the measured energy. The counts in each pixel were divided by the pixel area and the integration time and 
then interpolated onto a 500 X 500 image grid using as weights the inverse angular separation between image and data pixels out to a maximum separation equal 
to about 80% of the PSF. This image cube was summed over energy, and the color bar indicates the approximate rates in photons s sr~' with a square-root 
scale. Left: Photons selected by 0.08 < 4> < 0.20 and 0.56 < <f> < 0.80 for PSR J2032+4127 (pulsar "on"). Middle: Photons selected by 0.20 < 4> < 0.56 and 
0.80 < 4> < 1.08 (pulsar "off"), showing the background emission, including the two bright pulsars J2021+4026 and J2021+3651. Right: On-off difference 
map. The only remaining source is PSR J2032+4127. The 5° X 5° inset shows the PSF for a source with the same spectral energy dependence as measured for 
PSR J2032+4127, indicating that the pulsar is consistent with a point source. 



TABLE 1 

Source Positions Coincident with Be Star MT91 213 



Source Instrument R.A. (J2000.0) Decl. (J2000.0) Reference 



MT91213 MDM 2.4 m RETROCAM 20 32 13.137(18) +41 27 24.28(20) This work 
X-ray Chandra ACIS-I 20 32 13.143(24) +41 27 24.54(27) This work 
Radio GMRT 610 MHz 20 32 13.092(39) +41 27 24.16(48) Marti et al. (2007) 
PSR J2032+4127 Fermi LAT/GBT timing 20 32 13.07(4) +41 27 23.4(2) This work 

NOTE. — Units of right ascension are hours, minutes, and seconds, and units of declination are degrees, arcminutes, and arcseconds. 



w 1" FWHM. We have examined the area of the original ra- 
dio gridding localization of PSR J2032+4127 (§ l2~TlT > and 
find several X-ray sources that are each associated with an 
optical counterpart; the brightest is Cyg OB2 4, an R = 10.2 
07 III((f)) star (see Figure |7). Furthermore, this is a re- 
gion of enhanced diffu se hard X-ray emission identified by 
iMukheriee et al.l d2007l) . Figure [8] displays this ~ 1' feature 
and shows the timing position of PSR J2032+4127 for refer- 
ence. 

This area was surveyed for optical counterparts using the 
M DM Observatory Hiltne r 2.4 m telescope on 2002 August 
23 ( IMukheriee et al.ll2003l) . We compared the X-ray sources 
and their optical counterparts in Figure|7]to register the X-ray 
image with a correction of -0"16 in R.A. and +0."31 in decl. 
for the X-ray sources. Thus, we confirm that an X-ray source 
coincides with the radio position of PSR J2032+4127/GMRT 
5 and the star MT91 213 to within < 0."6, as listed in TableQ] 

MT91 indicate a spectral type B0 Vp and V = 1 1 .95 for star 
213. We obtained a spectrum of MT91 213 on 2009 June 13 
using the MDM 2.4 m telescope with Modspec, covering the 
wavelength range 4220-7550 A at 4 A resolution. Figure [9] 
shows the red portion of the spectrum, where emission lines 
of Ha and He I are evident. The equivalent width of Ha is 
-12.6 A, and the He I lines have asymmetric peaks that are 
separated by w 500 km s -1 , properties that are typical of Be 



stars. 

Since their spectra covered only the blue region, where 
emission lines are not obvious, MT9 1 did not report a Be clas- 
sification for star 213. Instead, they remarked that the absorp- 
tion lines appear unusually broad, and concluded that the star 
is multiple, thus the peculiar classification. Since all of the 
absorption lines in our spectrum are filled in with emission, 
we have no additional evidence to support the sug gestion that 
MT91 213 is a multiple system. In addition, Kimin kTet alj 
d2007l) did not find any evidence for binary motion among 
10 spectra of MT91 213 ob tained between 1999 and 2004. 
Massev & Thompson! d 19911) derive a color excess E(B-V) = 
1.43 for star 213, corresponding to Ay = 4.28. This opti- 
cal extinction would translate to X-ray co lumn density Nq = 
7.7 x 10 21 cm -2 using the conversion of Prede hT& SchmitJ 
(1995). We performed spectral fits to the 83 X-ray pho- 
tons detected from this source by Chandra. Power-law and 
Raymond-Smith plasma models fit equally well, with photon 
index T = 2.1 ±0.7 and kT = 4^2keV, respectively. The fit- 
ted A/h is consistent with the optical extinction to MT91 213. 
The unabsorbed X-ray flux in the 0.5-10keV band is Fx ~ 
3.2 x 10~ 14 erg cm -2 s -1 for the power law. The corresponding 
lumi nosity is Lx ~ 1.1 x 10 31 af 2 7 erg s -1 . 

Pare des et alj d2007l) suggested that the radio and X-ray 
sources coincident with MT91 213 are produced in a collid- 
ing wind binary. While we now know that the radio source 
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FIG. 7. — Left: Portion of Chandra image centered on the PSR J2032+4127 timing position, shown as a 3<r error ellipse. The Chandra 0.3-8 keV ACIS-I 
image has been exposure-corrected and smoothed with a 0"5 Gaussian kernel to highlight point sources. All X-ray point sources have an identified optical 
counterpart, shown on the right panel, including CXOU J203213. 5+41271 1, which is located 3"7 to the west of the bright star Cyg OB2 4. The X-ray source 
coincident with the star MT91 213 (see TablefT} may be coming either from the star or from PSR J2032+4127. Right: An S-band image of the X-ray field, taken 
with the MDM Observatory Hiltner 2.4 m telescope. The combined exposure time is 20 minutes, and the seeing is 1"0 (see Mukherjee et al. 2003 for further 
details). Cyg OB2 4, which masks the location of CXOU J203213. 5+41271 1, MT91 213, and other bright stars are saturated. Locations of Ch andra sources, 
astrometrically corrected using this image, are marked. The inset shows another image that was specially obtained in seeing of Of 6 (see § 12.2.41 to identify the 
optical counterpart of CXOU J203213.5+41271 1, 3"7 to the west of Cyg OB2 4. 
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FIG. 8.— Diffuse X-ray emission in the vicinity of PSR J2032+4127. The 
point sources have been removed from the Chandra 0.3-8 keV ACIS-I image 
of Figure[7J which is exposure-corrected and smoothed with a Gaussian ker- 
nel of <t — 14" to highlight diffuse emission. This image is a 7' X 7' zoom-in 
on Figure 3 from Mukherjee et al. (2007). The cross is located at the pulsar 
timing position (see Figure|7J. The dot-dashed circle approximates the extent 
of the HEGRA source TeV J2032+4130. 



is a pulsar, it is not clear what the relationship is between the 
radio pulsar, the X-ray source, and the Be star. The pulsar 
lacks any evidence of binary motion in its timing. The kine- 
matic contribution to the pulsar's frequency derivatives from 
its acceleration in orbit around a M c ~ 20M Q companion (ne- 



From Table [2] we obtain /// = -1.4 x 10" 13 s" 1 and /// = 
-2.4 x 10~ 22 s" 2 . If the pulsar is a binary companion of the Be 
star, the orbital period is probably P m ^ > lOOyr in order not 
to exceed these observed limits on the pulsar's acceleration. 

Such long-period systems are not known to exist (the 
longest is 5 years), presumably because such a wide binary 
would not survive the supernova kick. Instead, the pulsar is 
probably superposed by chance on the position of the star in 
this crowded field, in which case the X-ray source may be- 
long either to the radio pulsar or to the Be star, or both. The 
degraded PSF at the off-axis location of the Chandra image 
does not allow us to resolve this question. But the X-ray lu- 
minosity of the source is a fraction ~4x 10~ 5 £/ 2 7 of the pul- 
sar's E, which is wit hin the range of ob served efficiencies of 
young pulsars (e.g., Camilo et al. 2006). It is also compatible 
with the X-ray l uminosities of single early Be stars, although 
on the high end dBerghofer et al.lll997tlCohen et al.lfT997h . 

For completeness, we note that the location of one X- 
ray source that lies 3 "7 to the west of Cyg OB2 4, 
CXOU J2032 13. 5+4 12711, was saturated in all of our op- 
tical images, as shown in the right main panel of Figure [7] 
Its corrected coordinates are R.A. = 20 h 32 m 13?50, decl. = 
+41°27'11."8, with an error radius of w 0"3. We investigated 
the nature of this source because it is the only other hard X- 
ray source (with half of the 55 detected photons above 2 keV) 
close to the pulsar timing position. 
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FIG. 9. — Portion of the optical spectrum of MT91 213, showing emission lines that are typical of Be stars. The equivalent width of Ha is - 
He I lines have asymmetric peaks that are separated by 500kms~'. Interstellar absorption lines are marked below the spectrum. 



12.6 A, and the 



To identify CXOU J203213.5+41271 1 it was necessary to 
obtain an unsaturated optical image in good seeing. Images of 
1 s exposure were acquired on the MD M 2.4 m telesc ope on 
2009 May 24 using RETROCAM dMorganet al J 12001 and 
SDSS r filter in seeing of 0"6. The inset of Figure [7] (right) 
is the sum of 12 such exposures, which clearly reveals a star 
to the west of Cyg OB2 4 that we measure to fall within 0"3 
of the position of CXOU J203213. 5+41 271 1. On 2009 May 
25, we obtained an optical spectrum of this star on the MDM 
2.4 m telescope, using a 1" wide slit that cleanly isolated it. 
The resulting spectrum is that of a Be star, with an Ha emis- 
sion line as well as H and He absorption. 

The optical spectra of CXOU J203213.5+41271 1 and 
Cyg OB2 4 obtained on the same night show interstellar ab- 
sorption features that can be used to assess their relative dis- 
tances. The equivalent widths of the features agree to within 
~ 15%, comparable to their uncertainties, which implies 
that CXOU J203213.5+412711 is probably located in the 
Cyg OB2 association. There are too few X-ray counts from 
CXOU J203213. 5+412711 to obtain a spectrum, but we can 
estimate its flux by assumin g an absorbed power-law model. 
Massey & Thompson! d 19911) derive a color excess E(B— V) = 
1.59 for Cyg OB2 4 corresponding to A v = 4.50. This 
optical extinction would translate to X-r ay column density 
Nu = 8 x 10 21 cm" 2 using the conversion o flPredehl & Schrnita 
(1995). Assuming this JV H and T = 2.0, the unabsorbed 
flux of CXOU J203213.5+41271 1 in the 0.5-10keV band is 
Fx w 2.5 x 10~ 14 ergcm" 2 s" 1 . The corresponding luminosity 
is Lx « 9 x 1 30 d\ 7 erg s" 1 . X-ray luminositi es of Be stars de- 
cline sharply from early t o late spectral type (Berghof er et al.l 
1 1997b ICohen et all 1 19971) . The optical faintness of this star 
would seem to indicate a late spectral type. More precise 
spectral classification of CXOU J203213.5+41271 1 is needed 
to evaluate whether its high X-ray luminosity is anomalous. 
In any case, it is the second Be star discovered in the vicinity 
of PSR J2032+4127 that is apparently not connected to it. Cu- 
riously, this star and Cyg OB2 4 are both closer to the centroid 
of the diffuse emission in Figure [8] than is PSR J2032+4127, 
even though that may not have any physical significance. 

3. DISCUSSION 

With the detection of radio pulsations from two of the 
16 pulsars recently disco vered in blind sea rches of Fermi- 
LAT gamma-ray photons (lAbdo et al.ll2009fl) . we have taken 
a further step toward inferring the fraction of radio-quiet 
gamma-ray-emitting pulsars. Our non-detection of radio pul- 



sations from eight other LAT pulsars (§|2]l is not constraining 
enough, owing to insufficient sensitivity: minimum detectable 
flu x densities were S\a ~ 0.06 mJy for the three observed 
by iRoberts et al.l (120021) and S\a ~ 0.2 mJy for the remain- 
ing five, while very young radio pulsars have luminosities 
at least as low as L\ .4 = S\Ad 2 ~ 0.5mJykpc 2 (Cam ilo et al.l 
2002b), and the sma llest measured p ulsar luminosity is L\ 4 « 
0.03mJykpc 2 fsee lYoung etaD 19991). 

The Fermi LAT bright source list dAbdo etaii r2009c) cata- 
logs the 205 most significant sources above 0.1 GeV. Due to 
the spatially varying Galactic background, this is not a flux- 
limited sample, but we consider it here as a crude proxy for an 
unbiased sample containing a substantial number of pulsars. 

A total of 30 known pulsars are represented in the list, in- 
cluding 15 discovered blindly in LAT data. We exclude from 
consideration here three (previously known) millisecond pul- 
sars, since millisecond pulsars could have different emission 
properties. Of the remaining 27 pulsars, 13 are already known 
radio emitters, including the two reported here. Four oth- 
ers can be c onsidered to be rad i o quiet based on previous 
searches (see iBecker et al .ll2004t iHa toern et al.ll2004l 120071: 
iKassim & Lazioll 19991) . The status of the remaining 10 is still 
unclear, but given our success with PSRs J1741-2054 and 
J2032+4127, it is reasonable to suppose that some will even- 
tually be detected in the radio. Therefore, amo ng presently 
identif ied pulsars in the LAT bright source list dAbdo et al.l 
2009c), > 50% emit radio pulsations. 

In addition, however, nearly 50 of the 205 bright LAT 
sources are good pulsar candidates (they are without known 
AGN associations and are not variable). Regardless of 
whether gamma-ray pulsations can be detected from these, 
most of them are likely to be gamma-ray pulsars. Thus, in 
order to determine the true fraction of radio-quiet gamma-ray 
pulsars, it is crucial to do sufficiently sensitive radio searches 
of most of these 50 unidentified gamma-ray sources. 

We now discuss PSRs J1741-2054 and J2032+4127 in turn, 
for which we have gathered many of the measured and derived 
parameters in Table [2] 

3.1. PSRJ1741-2054 

The average radio flux density of PSR J 174 1-2054 is uncer- 
tain due to scintillation, but for our one detection at 1 .4 GHz, 
the nominal luminosity is L\ 4 = 0.025 li^mJykpc 2 , which is 
the smallest of any detected radio pulsar. Scintillation, which 
makes it difficult to observe this pulsar, is a result of the very 
small DM and thus small distance. At ~ 5 times the distance 



TWO RADIO-EMITTING GAMMA-RAY PULSARS 



9 



TABLE 2 

Measured and Derived Parameters for PSRs J1741-2054 and J2032+4127 



Parameter PSR J1741-2054 PSR J2032+4127 

Right ascension, R.A. (J2000.0) a 17 h 41 m 57?27(8) 20 h 32 m 13?07(4) 

Declination, decl. (J2000.0) a -20°53'13(34)" +41°27'23?'4(2) 

Galactic longitude, I (deg.) 6.42 80.22 

Galactic latitude, b (deg.) +4.90 +1.03 

Spin period, P (s) 0.41369996385(6) 0.14324743146(2) 

Period derivative, P (I0r [4 ) 1.694(3) 2.0063(1) 

Epoch (MJD) 54826.0 54840.0 

Timing data span (MJD) 54647-55004 54647-55033 

Timing residual, rms (ms) 3.4 0.36 

Dispersion measure, DM (pccirr 3 ) 4.7(1) 114.8(1) 

Rotation measure, RM (radirr 2 ) •■■ 215 ± 1 

Flux density at 2 GHz, S 2 (mJy) ■ ■ ■ b 0.12 ±0.03 

Radio-gamma-ray profile offset, S (P) 0.29±0.02 0.15±0.01 

Gamma-ray profile peak-to-peak separation, A (P) 0.23 ± 0.02 c 0.50 ± 0.01 

Gamma-ray (> 0.1 GeV) photon index, T 1.4 ±0.1 ±0.1 1.1 ±0.2 ±0.2 

Gamma-ray cut-off energy, E c (GeV) 1 . 1 ± 0.2 ± 0.2 3.0 ± 0.6 ± 0.7 

Photon flux (> 0.1 GeV) (lO" 8 cm" 2 s-') 20±1±3 7±1±2 

Energy flux (> 0. 1 GeV), P 7 (10" 1 1 erg cm -2 s _1 ) . 12±1±2 9±1±2 

Spin-down luminosity, E (erg s" 1 ) 9.4 X 10 33 2.7 X 10 35 

Characteristic age, t c (Myr) 0.39 0.11 

Surface dipole magnetic field strength (10 12 Gauss) 2.7 1.7 

Distance, d (kpc) d 0.4 3J> 

NOTE. — Values in parentheses are nominal TEMPO uncertainties on the last digit for parame- 
ters determined in timing fits. For gamma-ray spectral parameters, the first uncertainty is statistical 
and the second accounts for systematics. 

a The decl. of PSR J1741-2054 is known with greater precision than from this timing fit from 
a Swift observation (see Figure |2). In the timing fit for PSR J2032+4127 we account for 
rotational instability in the pulsar, parameterized by the second derivative of spin frequency: 
/= -1.7(3) X 10~ 21 s~ 3 . b The flux density received f rom P SR J1741-2054 varies greatly at 
1.4 GHz and 2 GHz, due to interstellar scintillation (see § 12.1.11 . This is the separation between 
PI and P3 (see Fig ure[3). d Distances are estimated from the DM and the electron density model 
of Co rdes & Lazio (2002). PSR J2032+4127 may be located at approximately half this distance 
(see §|3.2t. 

(still at only ~ 2 kpc), a pulsar of such luminosity would not parameters depend principally on the viewing angle ( (mea- 
be detectable with any existing telescope. This is worth con- sured from the rotation axis of the pulsar). For PSR J1741- 
sidering when discussing "radio quiet" neutron stars. 2054, the measured (A, 5) indicate a smaller £ than for most 

There are only eight pulsar s known with DM sma ller than other known gamma-ray pulsars. In detail, however, basic 
that of PSR J1741-2054 (see Manchester et al. 2005) 24 . Two geometric models still cannot explain many of the observed 
of them are exceedingly bright, and were discovered more (A, $), including for PSR J1741-2054, let alone features like 
than 40 years ago. Four others are millisecond pulsars, for the three peaks identified in this pulsar (Figure[3]and § 12.1.21 . 
which such a small DM is relatively easy to distinguish from and f urther ingredients are likely needed (see Watte rset al.1 
zero in searches, and thus from terrestrial radio interference. 2009). 

Only the remaining two were discovere d in modern searches In order to compute the flux of a gamma-ray pulsar aver- 

dTauris et al.1 fl994l lYoung etaTI [T999). Thus, while very aged over the sky, we need to know the geometry-dependent 
nearby pulsars are intrinsically rare, it is also the case that "beaming" factor /q that corrects the phase-averaged flux 
they have proven difficult to detect in radio surveys, with im- measured by Fermi. It is still not possible to infer this pre- 
plications for the completeness of the known sample. cisely for a given pulsar, especially when we do not pos- 

PSR J1741-2054 was d etected by EGRET as 3EG J1741- sess independent geometric information (such as from ra- 
2050 ( lHartman et aljfl999l) . Its EGRET flux above 0. 1 GeV dio pulsar polarization measurements or from high-resolution 
is similar to the LAT flux (§ 12.1.31 ). With P = 413 ms, r c = X-ray observations of pulsar wind nebulae). Nevertheless, 
0.39 Myr, and d w 0.4 kpc, PSR J 1741-2054 is reminiscent of most promising models — which i n addition to OG incl ude 
the nearest middle-aged pulsars, Geminga (P = 237 ms; t c = "two-pole caustic" (TPC) models ( iDyks & R udak 20QH) — 
0.34 Myr; d = 0.25 kpc) and PSR B0656+ 14 (P = 384 ms; t c = predict large fa. For the observed profile characteristics 
0.1 1 Myr; d = 0.29 kpc), although it has only 1/4 of their spin- of PSR J 174 1-2 054, which resemble somewhat those of 
down luminosity. " ' PSR B 1706^14 dPellizzoni et all l2009t) . the OG and TPC 

Among the known population of young (non-millisecond) mode ls appear to suggest that /q ~ 0.5-1 dWatters et al.l 
pulsars, PSR J1741-2054 has one of the smallest gamma- |2009|) . 

ray peak separations (A = 0.23) and largest gamma-ray-radio Based on the Fermi spectral parameters reported in § 12.1.31 

lags (5 = 0.29). Together, these values are roughly in agree- for PSR J1741-2054, we obtain L 7 (> 0. 1 GeV) = A-KfnF^d 2 = 
ment with the expectation from "outer-gap" (OG) gamma-ray 2.4 x I0 3i f a d^ 4 ergs' 1 =0.25 fa dl 4 E. Plausibly fa w 0.5- 
beam models, which predict an inver se dependence of A on i ; w hich might eventually be improved with further model- 
5 (e.g., | Romani & Yadigaroglu |[T995|). In OG models these ing . and ^ „ L PSR J1741-2054, which has the smallest E 

of any established non-millisecond gamma-ray pulsar, has a 

-■* http://www.atnf.csiro.au/research/pulsar/psrcat 
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high inferred efficiency for converting spin-down luminosity 
into gamma-rays, rj=Lj/E = 0.25 fn df )4 . This is in line with 
the tr end that rj increases with decreasing E (Thompso n et al.l 
1999). Constr aining the distance with VLBI measurements 
will be difficult, owing to the very small flux, but the proper 
motion could be accessible via either VLBI or HST mea- 
surements. Sensitive X-ray observations will also prove very 
useful in constraining the properties of this neutron star (see 

sua. 

Non-millisecond gamma-ray pulsars typically have ob- 
served pulse profiles with two distinct peaks. In many such 
cases, the trailing peak has a harder spectrum than PI. Al- 
though PSR J1741-2054 appears to have three distinct peaks, 
it follows this t rend i n that its trailing peak has a harder 
spectrum (see § 12.1.31 1. One respect in which PSR J 1 74 1 — 
2054 appears to stand out from non-millisecond pulsars is in 
having the smallest spectral cut-off energy reported so far, 
E c = 1.1 GeV. It remains to be determined whether this por- 
tends an incipient trend. 

3.2. PSR J2032+4127 

The radio luminosity of PSR J2032+4127, scaled to the 
standar d freq uency of 1 .4 GHz using its estimated spectral in- 
dex (§ EEQ, is L1.4 = 3dj 6 mJykpc 2 . While ~ 100 times 
larger than that of PSR J1741-2054, this is still a small lumi- 
nosity among the ob served young pulsar population (see, e.g., 
Camilo et al. 2002a). The actual distance to the pulsar may 
differ significantly fro m the 3.6 kpc estimated from the DM 
(Cord es & Lazioll2002l) . As we argue below, it may be that 

flf 3 . 6 ~ 0.5. 

As shown in Figure [4] we have made polarimetric observa- 
tions of PSR J2032+4127, and have measured Faraday rota- 
tion amounting to RM = +215radm~ 2 . This implies an aver- 
age Galactic magnetic field along the line of sight, weighted 
by the free electron density, of 2.3 /iG, which is a typical 
Galactic value. There are very few pulsars with measured 
RM in the approximate direction of this pulsar (lHan et al.1 
120061: lAbdo et al . 2009a), so that this measurement cannot be 
put into context and it does not provide a constraint on the 
pulsar distance. Also, we have tried to obtain information 
on the pulsar geometry by fitti ng a "rotating vector model" 
(Radhakris hnan & Cookd 1969) to the swing with pulse phase 
of the position angle of linear polarization (PA; top sub-panels 
of Figure 0). Unfortunately the fits are unconstrained, be- 
cause of the limited longitude coverage and relatively shal- 
low PA swing. The radio profile of PSR J2032+4127 ap- 
pears similar to those of several young pulsars in being sim- 
ple, relatively narrow, highly li nearly polarized, and in hav- 
ing little variation in PA (s ee Hohnston & Weisberg 120061: 
IWeltevrede & Johnstonll2008l) . 

The gamma-ray profile of PSR J2032+4127 (Figure [5} is 
broadly similar to t hose of several young pulsars, such as 
PSRs J2021+3651 dAbdo et al J l2009ah lHaipern et alj|2008t) 
and J0205+6449 (lAbdo et alJ l2009bl) . It has two narrow 
gamma-ray peaks preceded by the radio pulse, with (A,S) = 
(0.50,0.15). These values can be approximately understood 
within the context of both the OG and the TPC beam mod- 
els. In detail, however, S is somewhat larger than expected 
for t he observed A, giv en the simplest geometric models 
(see IWatters et al.l H009). This interpretive problem arises 
as well for PSRs J2021+3651 and J0205+6449. Unlike for 
PSR J1741-2054, these (A,<5) indicate a large viewing angle 
C for PSR J2032+4127, and a broad fan-like beam with cor- 



rection factor /q ~ 1. 

Based on the spectrum of PSR J2032+4127 (§ |2~23T >, we 
obtain L 7 (> 0.1 GeV) w 1.4 X lO^/n^.gergs" 1 . With/h ~ 1 
as indicated above, r/ w 0.5 d^ 6 . This is a very large nominal 
efficiency for a pulsar with E = 2.7 x lO ^ergs" 1 , whi ch has 
an open field line voltage of 10 I5 V (cf. lAronsl[l996h . Part 
of the answer may lie in a sm aller d istance, and indeed it is 
plausible that c/3.6 ~ 0.5 (see § 13.2. 11 1. 

PSR J2032+4127 is located in the Cygnus region, which 
contains several bright gamma-ray point sources and strong 
spatially varying diffuse emission. This leads to particular 
difficulty in modeling the spectrum of the pulsar (see § 12.2.31 
and Figure [6). PSR J2032+4127 is marginally consistent 
with the location of the EGRET source 3 EG J2033+4118 
dHartman et al.ll 19991) . which is listed as EGR J2033+41 17 i n 
the revised EGRET catalog of Casandj ian & Grenierl (120081) . 
We therefore consider PSR J2032+4127 as the likely iden- 
tification of the EGRET source. The photon fluxes above 
0.1 GeV for 3EG J2033+41 18 and EGR J2033+41 17 are, re- 
spectively, 73 ± 7 and 5 2 ± 7, in units of 1 0~ 8 cm -2 s" 1 . In the 
Fermi bright source list (lAbdo et al.l2009d) . PSR J2032+4127 
corresponds to source 0FGL J2032. 2+4122, which has a flux 
of 54 ± 5. In § 12331 however, we report for PSR J2032+4127 
a flux of (7 ± 1 ± 2) x 10" 8 cm" 2 s" 1 , or a bout 1/8 of t h at give n 
for 0FGL J2032.2+4122. Compared to lAbdo et all d2009d) . 
we have used significantly more data, an updated LAT re- 
sponse function, an improved diffuse ba ckgrou nd model , con- 
sidered only on-pulse emission (see §§ 12. 2. 31 and 12.1.31 ). and 
modeled additiona l point sources in the field. Finally, the 
lAbdo et all (l2009d) spectral model was a simple power law, 
generally inadequate to describe the spectral breaks found in 
pulsars. Together, these differences presumably account for 
the large discrepancy in photon flux. The new, much smaller, 
flux value nevertheless implies an unreasonably large effi- 
ciency T] w 0.5 d\ 6 . We believe that for PSR J2032+4127 our 
flux estimate is the most reliable one available at this point, 
but still subject to improvement as the modeling of the com- 
plex Cygnus region evolves. And we also suggest that a rea- 
sonable conversion efficiency for PSR J2032+4127 is most 
naturally achieved with dj (, w 0.5. 

The traili ng pea k of PSR J2032+4127 has a harder spectrum 
than PI (§ 12.2.31 ). This behavior is si milar to that observ ed 
in PSRs J202 1+3651 and J0205+6449 (lAbdo et al.ll2009allbh . 
which as already noted have pulse profiles broadly similar to 
PSR J2032+4127, but have E that are larger by factors of 10 
and 100, respectively. This gathering trend now applies to 
many known gamma-ray pulsars with di stinct peaks from 
PSR J1741-2054 (§ IXOb to the Crab (iThompsonl l200l . 
spanning a factor of 5000 in E. 

3.2.1. The (formerly) unidentified TeV J2032+4130 

TeV J2032+4130 was discovered serendipitously by the 
HEGRA system of Atmospheric Cherenkov Telescopes in 
obser vations of the Cygnus region over the period 1999- 
2001 dAharonian et all [2002). It was the first unclassified 
TeV source, and its origin remained undetermined until now. 
Analysis of combined HEGRA data from 1999-2002 gave 
a final position centroid for the extended TeV source of 
R.A. = 20 h 31 m 57?0±6. s 2 sm ,± 13?7 sys , deck = +41°29'57" ± 
1 fl s tat ± UOsvs, and a Gaussia n radius of a = 6. '2 ± 1 f2 stat ± 
0.'9 sys dAharonian et al.ll2005h . The HEGRA-measured pho- 
ton flux above ITeV is (6.9 ± 1.8) x 10" 13 cm -2 s" 1 with 
T = 1.9 ±0.3. MAGIC detected TeV J2032+4130 at a 
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best position of R.A. = 20 h 32 m 20 s ± 11^, ± 11^, decl. = 
+41°3 0'36" ± 1 .'2 stat ± 1 f 8 sys , and radius a = 6.'0 ± l.'7 stat ± 
0f6 sys (I Albert et al.l |2008|). The flux measured by MAGIC 
above 1 TeV is 4.5 x 10~ 13 cm" 2 s _1 with photon index V = 
2.0 ± 0.3. The Whipple Observatory detected a source of lu- 
minos ity = 4 x lO 33 ^ 2 -, erg s" 1 , assuming a Crab-like spec- 
trum ( Konopelko et al. 2007), which is about a factor of 2 
greater than measured by HEGRA, but probably consistent 
given the uncertainty in the spectrum. The Whipple position 
and extent of the source are marginally consistent with those 
of HEGRA. Milagro detects a diffuse source in a 3° x 3° re- 
gion centered on Te V J2032+4130, th at exceeds the HEGRA 
flux by a factor of 3 (lAbdo et al.ll2007h . This region may con- 
tain multiple source s; the Milagro ex cess at the location of 
the pulsar is 7.6 a dAbdo et al.ll2009el) and is consistent with 
an extrapolation of the HEGRA spectrum. None of the ex- 
periments have found strong evidence for flux variability of 
TeV J2032+4 1 30 in the period 1 999-2005 . 

Pulsar wind nebulae (PWNe) comprise the largest class of 
Galactic TeV sources 25 , and all but the youngest, with ages 
of ~ 10 3 yr, are spatially extended. PSR J2032+4127 is lo- 
cated within the 1 a extent of TeV J2032+4130 (see Figure©, 
only 4' from its HEGRA centroid. We therefore propose its 
PWN as the source of TeV J2032+4130. Further support for 
the association of PSR J2032+4127 with TeV J2032+4130 
comes from comparing their properties with those of other 
TeV PWNe. 

First, however, we discuss the implications of the appar- 
ent coincidence of PSR J2032+4127 with the massive young 
star cluster Cyg OB2. Distance estimates to Cyg OB2 
range over 1.45-1.7kpc (lHansonl20 03; Mass ev & Tho mpson 
11991b . while the nominal DM distance of PSR J2032+4127 is 
3.6 kpc. Considering the possibility of unmodeled local en- 
hancements in the electron density that may cause the DM 
model to overestimate the pulsar distance, we do not rule 
out that PSR J2032+4127 could be colocated with Cyg OB2 
(which, as noted in § 13.21 would imply a far more reason- 
able pulsar GeV conversion efficiency r\). PSR J2032+4127 
is projected 15' from the center of the cluster, which has a 
half-light radius of 13' an d a diameter of ~ 2° as determined 
from 2MASS star counts (lKnodlsededl2000l) . 

Conventional wisdom holds that no neutron stars have been 
born in Cyg OB2 because the age of the majority of its 
stars, determined from isochrone fittin g, is only 2-2.5 Myr 
dHansonll20"03l: iNegueruela et al.ll2008l) . At this age, stars of 
M < 35 Mq are still close to the main sequence, and no super- 
novae have occurred. However, there is evidence for earlier 
episodes of star formation in the immediate nei ghborhood, 
such a s the ~ 5-7 Myr-old A stars discovered by iDrew et all 
(2008) within a 1° radius of Cyg OB2 and at the same dis- 
tance. Ol der OB s t ars co incident with Cyg OB2 are also dis- 
cussed by Hanson (2003). Therefore, we cannot rule out that 
the progenitor star of PSR J2032+4127 was born in a recent 
episode of star formation at the same distance as Cyg OB2, 
while the typical neutron star velocity of ~ 250km s _I would 
allow it to have travelled ~ 1° at that distance in 10 5 yr. 

The offset of TeV J2032+4 1 30 from the center of Cyg OB2 
is most easily explained by the existence and location of 
PSR J2032+4127, i.e., that PSR J2032+4127 is the source 
of TeV J2032+4130, and not, e.g., winds from O stars in 
the cluster. In the latter case, it would be difficult to under- 

25 http://www.mppmu.mpg.de/~rwagner/sources 



stand why the TeV emission is not more widely distributed 
among the dozens of O stars of Cyg OB2, and centered on 
the cl uster, which is quite spherical in 2MASS (Kntidlseder 
120001) . This argument for associating PSR J2032+4127 with 
TeV J2032+4130 applies whether or not PSR J2032+4127 is 
actually in the cluster. 

In order to compare the efficiency of TeV gam ma-ray pro- 
ductio n with those of other pulsars listed by iGallant et al.l 
(2008), we integrate the HEGRA spectrum over 0.3-30 TeV, 
giving F(0.3-30 TeV) = 5.1 x lO-^ergcnrV 1 , L 7 = 1.8 x 
10 33 flf 2 7 ergs -1 , and e = L 1 /E = 0.007t/ 2 7 . This is consis- 
tent with the range of efficien cies e = IP -4 — 0. 1 1 found for 
other PWNe and ca ndidates by IGallant et all (120081) (see also 
iHessels et al.ll2008l) . Even at d = 3.6 kpc, e = 0.03 is not ex- 
ceptional. Also typical for a PWN is the flux of the diffuse X- 
rays within 1' of PSR J2032+4127, ~ 1 x IP" 13 erg cm" 2 s" 1 
in the 0.5-10keV band (see Figure^ Mukherjee et al.ll2007l) . 
which corresponds to Lx/E ~ 1 x 10~ 4 c/ 2 7 . Therefore, we 
approximate the ratio L-y/Lx ~ 50, which is compa rable to 
the ratio for other older pulsars (Matta na et al.ll20()9). How- 
ever, we do not see any evidence for a larger X-ray nebula 
co vering almos t the e ntire HEGRA source that was claimed 
bv iHornset ail (120071) using XMM images, with an order of 
magnitude higher flux than our small nebula. 

The spin-down luminosity of PSR J2032+4127 stands out 
as the smallest of any pulsar identified with a TeV source; 
the others have E > 10 36 ergs _I . One possible exception is 
PSR J 1702-4 128 with E = 3.4 x 10 35 erg s" 1 , if it is associ- 
ated with HESS J1 702-420. The latter id entification is prob- 
lematic, however ( Aharo riian et al.l 120081) . because it would 
require e ~ 0.7 at the estimated distance of 5 kpc, and a rather 
extreme asymmetry for the TeV emission, peaking more than 
0.°6 from the pulsar. Nevertheless, the existence of still 
unidentified TeV sources allows the possibility that new pul- 
sar identifications may be found that are less energetic than 
the already known pulsar counterparts, and TeV J2032+4130 
is probably one such example. 

Sho rtly after the discovery of TeV J2032+4130. lBednarekl 
(2003) proposed that 3EG J2033+4118 is a Vela-type pul- 
sar in Cyg OB2, and that TeV J2032+4130 is its PWN. 
Using a time-dependent model for the PWN that includes 
both hadronic and leptonic processes, he assumed a birth 
period of 2ms, 5~6x 10 12 G, and present pulsar param- 
eters P = 210 ms and t c = 2 x 10 4 yr, corresponding to E = 
2.3 x 10 35 ergs _1 . The latter is close to the measured value 
for PSR J2032+4127 of E = 2.7 x 10 35 ergs" 1 , although the 
measured characteristic age of 1.1 x 10 s yr suggests that the 
time-dependence of the model is not well constrained. For 
example, PSR J2032+4127 may be too old for its initial spin 
energy to h ave had much influence on its p resent TeV lumi- 
nosity. See de Jage r & Diannati-Atail ((2008) for a discussion 
o f relevant timesca les . 

iButt et ail d2008l) point to a possible shell-like arrangement 
of predominantly non-thermal radio emi ssion, centered on th e 
TeV source, in VLA images (see also Par edes et al.l l2007h . 
While it is not even clear that this is a single coherent struc- 
ture, it is unlikely to be the remnant of the supernova that gave 
birth to the ~ 10 5 yr-old PSR J2032+4127, as its radius of 3 pc 
at the distance of Cy g OB2 would be to o small. Using higher 
resolution VLA data. iMartf et al.l d2007l) do not confirm a su- 
pernova remnant nature for this emission. 

Although PSR J2032+4127 is consistent in position with 
the Be star MT91 213, there is no other evidence that they are 
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a binary pair such as the prototype gamma-ray emitting binary 
PSR B1259-63, which has an orbital period of 3.4 yr. Never- 
theless, if the proximity of PSR J2032+4127 to Cyg OB2 is 
real rather than apparent, it is likely that the local radiation 
background from the massive star cluster enhances inverse 
Compton TeV emission from the PWN, which might other- 
wise not have been detectable. Although the total luminos- 
ity from massive stars in Cyg OB2 is uncertain, estimates of 
~ 100 O stars lead to a luminosity of ~ 10 41 ergs _1 , and a 
photon energy density of ~ 10 2 eVcm" 3 , which at the loca- 
tion of PSR J2032+4127 may manifest itself mainly as repro- 
cessed IR emission from dust (compared with 0.26 eV cm" 3 
for the cosmic microwave background). Detailed modeling, 
including inverse Compton scattering in the Klein-Nishina 
regime, will be needed to evaluate the possible TeV emis- 
sion using this enhanced background radiation, and could 
help constrain the distance between PSR J2032+4127 and 
Cyg OB2. 

4. CONCLUSIONS 

We report radio pulsar detections of two Fermi sources pre- 
viously known only as gamma-ray pulsars, as well as details 
of their gamma-ray spectra, and probable identifications of 
their X-ray counterparts. Gamma-ray efficiencies are esti- 
mated using their radio dispersion-measure distances. Both 
pulsars have hard power-law gamma-ray spectra with expo- 
nential cut-offs in the GeV band, and large efficiencies rela- 
tive to spin-down power of ~ 0.2, assuming isotropic emis- 
sion. The phase offsets between their radio and gamma-ray 
pulses follow trends observed in other pulsars, probably indi- 
cating their viewing geometry, and consistent with outer mag- 
netosphere models, although more physics must be added to 
existing calculations to model the detailed structure and spec- 
tral evolution of the observed light curves. 

PSR J1741-2054 has a low E = 9.4 x 10 33 ergs" 1 , and its 
characteristic age of 0.4 Myr is compatible with the soft X- 
ray spectral characteristics of its putative counterpart, inferred 
to be surface thermal emission plus a non-thermal compo- 
nent. With a radio luminosity smaller than that of any other 
known radio pulsar, PSR J1741-2054 at a distance of 0.4 kpc 
is approximately twice as far as the radio-quiet, middle-aged 
gamma-ray pulsar Geminga, which it resembles in its spin- 
down parameters and X-ray properties. PSR J1741-2054 be- 
gins to answer the question of where are the other Gemingas. 

PSR J2032+4127 is a more energetic pulsar, with E = 2.1 x 
10 35 ergs _1 , and is brighter in radio. A precise timing posi- 



tion derived from a joint fit to Fermi and GBT data confines 
the pulsar position to a region of diffuse X-ray emission previ- 
ously identified in a Chandra image, which is now presumed 
to be its PWN. The location of PSR J2032+4127 within the 
extent of the diffuse HEGRA source TeV J2032+4130 solves 
the 10-year-old mystery of the nature of this, the first unidenti- 
fied TeV source. PSR J2032+4127 is probably one of the least 
energetic pulsars powering TeV PWNe, which are now known 
to be the most numerous type of Galactic TeV source. The lo- 
cation of PSR J2032+4127 projected close to the core of the 
massive, young stellar association Cyg OB2 at a distance of 
1.5-1. 7 kpc suggests that this is its true distance, rather than 
3.6 kpc estimated from the radio pulsar dispersion. 



We thank profusely all the Fermi team members who have 
built and who operate the magnificent spacecraft and LAT ex- 
periment. The GBT is operated by the National Radio As- 
tronomy Observatory, a facility of the National Science Foun- 
dation operated under cooperative agreement by Associated 
Universities, Inc. The Parkes Observatory is part of the Aus- 
tralia Telescope, which is funded by the Commonwealth of 
Australia for operation as a National Facility managed by 
CSIRO. We are grateful to the Swift project scientist and staff 
for the observation of the LAT PSR J1741-2054 field, and to 
John Thorstensen for taking the MDM optical spectrum. 

The Fermi LAT Collaboration acknowledges generous on- 
going support from a number of agencies and institutes that 
have supported both the development and the operation of 
the LAT as well as scientific data analysis. These include the 
National Aeronautics and Space Administration and the De- 
partment of Energy in the United States, the Commissariat a 
FEnergie Atomique and the Centre National de la Recherche 
Scientifique / Institut National de Physique Nucleaire et de 
Physique des Particules in France, the Agenzia Spaziale Ital- 
iana and the Istituto Nazionale di Fisica Nucleare in Italy, the 
Ministry of Education, Culture, Sports, Science and Technol- 
ogy (MEXT), High Energy Accelerator Research Organiza- 
tion (KEK) and Japan Aerospace Exploration Agency (JAXA) 
in Japan, and the K. A. Wallenberg Foundation, the Swedish 
Research Council and the Swedish National Space Board in 
Sweden. 

Facilities: CXO (ACIS-I), Fermi (LAT), GBT (BCPM, 
GUPPI), Hiltner (RETROCAM), Parkes (PMDAQ), Swift 
(XRT) 



REFERENCES 



Abdo, A. A., et al. 2007, ApJ, 658, L33 

— . 2009a, ApJ, 700, 1059 

— . 2009b, ApJ, 699, L102 

— . 2009c, ApJS, 183,46 

— . 2009d, ApJ, 696, 1084 

— . 2009e, ApJ, 700, LI 27 

— . 2009f, Science, 325, 840 

— . 2009g, Science, 325, 848 

Aharonian, E, et al. 2002, A&A, 393, L37 

— . 2005, A&A, 431, 197 

— . 2008, A&A, 477, 353 

Albert, J., et al. 2008, ApJ, 675, L25 

Arons, J. 1996, A&AS, 120, 49 

Atwood, W. B., et al. 2009, ApJ, 697, 1071 

Backer, D. C, Dexter, M. R., Zepka, A., D., N., Wertheimer, D. J., Ray, P. S., 

& Foster, R. S. 1997, PASP, 109, 61 
Becker, W., et al. 2004, ApJ, 615, 897 
Bednarek, W. 2003, MNRAS, 345, 847 



Berghofer, T. W., Schmitt, J. H. M. M., Danner, R., & Cassinelli, J. P. 1997, 
A&A, 322, 167 

Butt, Y. M., Combi, J. A., Drake, J., Finley, J. P., Konopelko, A., Lister, M., 

Rodriguez, J., & Shepherd, D. 2008, MNRAS, 385, 1764 
Butt, Y. M., Drake, J., Benaglia, P., Combi, J. A., Dame, T., Miniati, E, & 

Romero, G. E. 2006, ApJ, 643, 238 
Camilo, E, Manchester, R. N., Gaensler, B. M., Lorimer, D. L., & Sarkissian, 

J. 2002a, ApJ, 567, L71 
Camilo, E, Ransom, S. M., Gaensler, B. M., Slane, P. O., Lorimer, D. R., 

Reynolds, J., Manchester, R. N., & Murray, S. S. 2006, ApJ, 637, 456 
Camilo, E, et al. 2002b, ApJ, 571, L41 
Casandjian, J.-M., & Grenier, I. A. 2008, A&A, 489, 849 
Cohen, D. H., Cassinelli, J. P., & MacFarlane , J. J. 1997, ApJ, 487, 867 
Cordes, J. M., & Lazio, T. J. W. 2002, arXiv:astro-ph/0207156 
de Jager, O. C, & Djannati-Atai, A. 2008, in Neutron Stars and Pulsars: 40 

Years After Their Discovery, ed. W. Becker (Berlin: Springer), 45 1 
Drew, J. E., Greimel, R., Irwin, M. J., & Sale, S. E. 2008, MNRAS, 386, 1761 
Dyks, J., & Rudak, B. 2003, ApJ, 598, 1201 



TWO RADIO-EMITTING GAMMA-RAY PULSARS 



13 



Gallant, Y. A., et al. 2008, in AIP Conf. Ser. 983, 40 Years of Pulsars: 
Millisecond Pulsars, Magnetars and More, ed. C. Bassa, Z. Wang, 
A. Cumming, & V. M. Kaspi (Berlin: Springer), 195 

Halpern, J. P., et al. 2008, ApJ, 688, L33 

Halpem, J. P., Camilo, R, & Gotthelf, E. V. 2007, ApJ, 668, 1 154 

Halpern, J. P., Gotthelf, E. V., Camilo, E, Helfand, D. J., & Ransom, S. M. 

2004, ApJ, 612, 398 
Han, J. L., Manchester, R. N., Lyne, A. G., Qiao, G. J., & van Straten, W. 

2006, ApJ, 642, 868 
Hanson, M. M. 2003, ApJ, 597, 957 

Harding, A. K., Grenier, I. A., & Gonthier, P. L. 2007, Ap&SS, 309, 221 
Hartman, R. C, et al. 1999, ApJS, 123, 79 
Hessels, J. W. T., et al. 2008, ApJ, 682, L41 

Horns, D., Hoffmann, A. I. D., Santangelo, A., Aharonian, F. A., & Rowell, 

G. P. 2007, A&A, 469, L17 
Hotan, A. W., van Straten, W., & Manchester, R. N. 2004, Proc. Astr. Soc. 

Aust., 21, 302 

Johnston, S., & Weisberg, J. M. 2006, MNRAS, 368, 1856 
Kassim, N. E., & Lazio, T. J. W. 1999, ApJ, 527, L101 
Kiminki, D. C., et al. 2007, ApJ, 664, 1102 
Knodlseder, J. 2000, A&A, 360, 539 
Konopelko, A., et al. 2007, ApJ, 658, 1062 

Lorimer, D. R., & Kramer, M. 2005, Handbook of Pulsar Astronomy 

(Cambridge University Press) 
Manchester, R. N, Hobbs, G. B., Teoh, A., & Hobbs, M. 2005, AJ, 129, 1993 
Manchester, R. N., et al. 2001, MNRAS, 328, 17 

Marti, J., Paredes, J. M., Ishwara Chandra, C. H., & Bosch-Ramon, V. 2007, 

A&A, 472, 557 
Massey, P., & Thompson, A. B. 1991, AJ, 101, 1408 (MT91) 



Mattana, E, et al. 2009, ApJ, 694, 12 

Morgan, C. W., Byard, P. L., DePoy, D. L., Derwent, M., Kochanek, C. S., 

Marshall, J. L., O'Brien, T. P., & Pogge, R. W. 2005, AJ, 129, 2504 
Mukherjee, R., Gotthelf, E. V., & Halpern, J. P. 2007, Ap&SS, 309, 29 
Mukherjee, R., Halpern, J. P., Gotthelf, E. V., Eracleous, M., & Mirabal, N. 

2003, ApJ, 589, 487 
Negueruela, I., Marco, A., Herrero, A., & Clark, J. S. 2008, A&A, 487, 575 
Paredes, J. M., Marti', J., Ishwara Chandra, C. H., & Bosch-Ramon, V. 2007, 

ApJ, 654, L135 
Pellizzoni, A., et al. 2009, ApJ, 691, 1618 
Predehl, P., & Schmitt, J. H. M. M. 1995, A&A, 293, 889 
Radhakrishnan, V., & Cooke, D. J. 1969, Astrophys. Lett., 3, 225 
Ransom, S. M. 2001, PhD thesis, Harvard Univ. 

Ransom, S. M., Eikenberry, S. S., & Middleditch, J. 2002, AJ, 124, 1788 
Roberts, M. S. E., Hessels, J. W. T., Ransom, S. M., Kaspi, V. M., Freire, 

P. C. C, Crawford, E, & Lorimer, D. R. 2002, ApJ, 577, L19 
Romani, R. W., & Yadigaroglu, I.-A. 1995, ApJ, 438, 314 
Tauris, T. M., et al. 1994, ApJ, 428, L53 

Thompson, D. J. 2001, in AIP Conf. Ser. 558, High Energy Gamma-Ray 
Astronomy, ed. F. A. Aharonian & H. J. Volk (Berlin: Springer), 103 

Thompson, D. J. 2004, in Cosmic Gamma-Ray Sources, ed. K. S. Cheng & 
G. E. Romero (Dordrecht: Kluwer), 149 

Thompson, D. J., et al. 1999, ApJ, 516, 297 

Watters, K. P., Romani, R. W., Weltevrede, P., & Johnston, S. 2009, ApJ, 695, 
1289 

Weltevrede, P., & Johnston, S. 2008, MNRAS, 391, 1210 

Young, M. D., Manchester, R. N., & Johnston, S. 1999, Nature, 400, 848 



